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Abstract 

Batch  and  miscible-displacement  experiments  were  conducted  to  determine  the 
extent  of  adsorption-desorption  and  transport  of  2, 4,6-trlnltrotoluene  (TNT)  and 
2, 3, 5-trinltro- 1 , 3, 5-trlazlne  (RDX)  In  soils,  A  reference  bentonite  clay,  contami¬ 
nated  (Kolin)  soil  from  the  Louisiana  Army  Ammunition  Plant,  and  two 
uncontaminated  soils  were  used.  The  TNT  Isotherm  for  bentonite  clay  was 
described  equally  well  using  linear,  Fieundllch,  Langmuir  and  modified  Langmuir 
models.  TNT  adsorption  and  desorption  Isotherms  showed  a  lack  of  hysteric 
behavior,  with  TNT  retention  as  a  fully  reversible  mechanism.  Transport  results 
from  bentonite  clay  columns  Indicated  that  TNT  was  a  highly  mobile  contami¬ 
nant  and  fully  conservative  In  the  presence  of  methanol  as  the  background 
solution.  Mobility  of  TNT  was  strongly  retarded,  with  some  50%  of  the  applied 
TNT  retained  In  the  bentonite  clay  column  when  0,005-M  Ca(N03)2  was  the 
background  solution.  Transport  results  In  Norwood  (fine  silty)  soli  columns 
revealed  that  TNT  was  strongly  retained  In  this  low-organic-matter  and  low-clay- 
content  soil.  The  use  of  a  transport  model  with  either  Freundllch  or  linear 
retention  and  an  irreversible  mechanism  predicted  the  TNT  transport  data  well, 
Poor  predictions  were  obtained  when  model  parameters  based  on  batch 
retention  data  were  used.  Reasons  for  model  failure  may  be  attributable  to  TNT 
sorptlon-klnetlcs  and  retention  because  of  diffusion  Into  clay  lattices,  Transport 
results  Indicated  high  mobility  of  RDX  with  limited  retardation,  which  Is 
consistent  with  release  and  transport  data  from  the  contaminated  soli, 


Cover:  Experimental  setup. 


For  conversion  of  SI  metric  units  to  U  ,S, /British  customary  units  of  measurement 
consult  ASTM  Standard  E380-89a,  Standard  Practice  for  Use  of  the  International 
System  of  Units t,  published  by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa,  19103, 


CRREL  Report  94-7 

U.S.  Army  Ccrps 
of  Engineers 

Cold  Regions  Research  & 
Engineering  Laboratory 

Sorption-Desorption  and  Transport  of 
TNT  and  RDX  in  Soils 


H.  Magdl  Selim  and  Iskandar  K,  Iskandar 


May  1 994 


Acceslon  For 

NTIS  CRA&I 
DTIC  TAB 
Unannounced  □ 
Justification 


By . 

Distribution/ 


Availability  Codes 


Dlst 


m 


Avail  and/or 
Special 


Prepored  for 

OFRCE  OF  THE  CHIEF  OF  ENGINEERS 


Approved  for  public  release;  distribution  Is  unlimited. 


PREFACE 


Tills  report  was  prepared  by  Dr.  H.  Magdi  Selim,  Professor  of  Soil  Physics  at 
Louisiana  State  University,  Baton  Rouge,  Louisiana,  and  Dr.  Iskandar  K.  Iskandar, 
Chief  of  Geochemical  Sciences  Branch,  Research  Division,  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory ,  This  work  was  carried  under  out  the  Strategic 
Environmental  Research  Development  Program  (SERDP),  Fate  and  Transport  in  Season¬ 
ally  Frozen  Soils  and  Discontinuous  Permafrost,  and  the  Department  of  Defense 
RDT&E  Research  Program, 

The  authors  thank  Janet  Ardoin,  former  Research  Associate,  for  her  technical 
assistance,  and  Paul  Miyares  and  Dr.  Thomas  F,  Jenkins,  both  of  CRREL,  for  their 
assistance  with  HPLC  methodology.  The  help  of  Susan  E,  Hardy  and  Arlene  Phillips 
in  editing  the  final  version  of  the  manuscript  is  appreciated. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional 
purposes.  Citation  of  brand  names  does  not  constitute  an  official  endorsement  or 
approval  of  the  use  of  such  commercial  products. 


ii 


CONTENTS 


Preface .  ii 

Introduction . 1 

Transport  models .  2 

Experimental,. . 4 

Soils  and  reagents . 4 

Sorption  experiments . 5 

Desorption  experiments .  5 

Transport  experiments .  5 

Results  and  discussion.... . . . 7 

Sorption  Isotherms . 7 

Desorption  and  hysteresis . 9 

Transport . 12 

Contaminated  soil . 20 

Summary  and  conclusions . 2:* 

Literature  cited. . 23 


ILLUSTRATIONS 


Figure 

1.  TNT  sorption  on  bentonite  clay . 

2.  TNT  sorption  on  bentonite  clay  and  fitted  Freundlich  and 

Langmuir  models . 

3.  TNT  sorption  on  bentonite  clay  and  fitted  modified  Langmuir 

model  ...in. . . . . . . . y . . . . 

4.  TNT  sorption-desorption  isotherms  for  bentonite  clay . 

5.  TNT  concentration  in  solution  in  the  bentonite  clay  vs.  number 

of  dilutions-extractlons  during  desorption . 

6.  Amount  of  TNT  sorbed  on  bentonite  clay  vs.  number  of 

dilutions-extractions  during  desorption . 

7.  TNT  sorption-desorption  isotherms  for  bentonite  clay  for  all 

initial  TNT  concentrations  used . 

8.  TNT  breakthrough  results  for  a  bentonite-sand  column  with 

input  TNT  concentration  of  100  Hg/mL  in  a  50%  methanol 
background  solution . 

9.  TNT  breakthrough  results  for  a  bentonite-sand  column  with 

an  input  TNT  concentration  of  10.65  ).ig/mL  in  a  0,005-M 
Ca(N03)2  background  solution . 

10.  RDX  breakthrough  results  for  a  bentonite-sand  column  with 

input  RDX  concentration  of  10.65  pg/mL  in  a  0.005-M 
Ca(N03)2  background  solution . 

11.  RDX  breakthrough  results  for  Kolin  uncontaminated  soil 

column  with  input  RDX  concentration  of  10.04  gg/mL 
in  0,005-M  Ca(NC)3)2 . 


8 

8 

9 

10 

11 

11 

12 

13 

14 

16 


17 


iii 


12.  TNT  breakthrough  results  for  a  Norwood  soil  column .  19 

13.  Desorption  or  release  of  RDX  and  TNT  from  the  contaminated  soil 

samples  obtained  from  the  Louisiana  Army  Ammunition  Plant .  21 

14.  TNT  and  RDX  breakthrough  results  for  the  contaminated  soil 

samples  obtained  from  the  Louisiana  Army  Ammunition  Plant .  22 


TABLES 

Table 

1.  Parameters  of  the  various  soil  columns  for  the  TNT  and  RDX 

miscible  displacement  experiments .  6 

2.  Concentration  of  TNT  and  RDX  input  pulses,  background  solution 

and  pulse  duration  for  the  miscible  displacement  column 

experiments . . .  6 

3.  Best-fit  model  parameters  for  batch  sorption  data  of  TNT  on  reference 

clay . . . 10 

4.  Best-fit  model  parameters  for  the  TNT  and  RDX  miscible  displacement 

experiments  . . 16 

5.  Concentrations  of  explosive  residues  released  from  contaminated  soil 

samples  from  the  Louisiana  Army  Ammunition  Plant .  20 


iv 


Sorption-Desorption  and  Transport  of 
TNT  and  RDX  in  Soils 

H,  MAGDI  SELIM  AND  ISKANDAR  K.  ISKANDAR 


INTRODUCTION 

When  the  Superfund  Program  was  enacted  in 
1980,  a  1.6-billion-dollar  trust  fund  was  autho¬ 
rized  for  the  cleanup  of  the  worst  hazardous  waste 
sites.  However,  in  1986,  the  Environmental  Protec¬ 
tion  Agency  (EPA)  estimated  the  need  for  five 
times  this  amount,  and  an  8.5-billion-dollar  fund 
was  approved.  In  addition,  the  Department  of 
Defense  (DOD)  estimated  24,5-billion  dollars  for 
cleanup  of  DOD  sites,  while  the  Department  of 
Energy  estimated  at  least  30-btllion  dollars  for 
these  sites  for  the  next  5  years.  EPA  (1993)  esti¬ 
mated  that  the  entire  cost  for  cleanup  will  exceed 
1-trillion  dollars. 

One  reason  for  this  uncertainty  In  the  cost  of 
cleanup  Is  the  lack  of  knowledge  about  the  envi¬ 
ronmental  fate  and  toxicity  of  pollutants  such  as 
organics,  heavy  metals,  Army-unique  compounds 
and  radioactive  nuclides,  and  their  effects  on  ter¬ 
restrial  and  aquatic  life  forms  as  well  as  on  hu¬ 
mans.  While  cleanup  is  required  by  regulations 
and  statutes,  there  are  no  et  i vironmental  criteria  or 
health  advisories  available  for  the  military  materi¬ 
als  of  interest.  The  ability  to  make  economical  and 
timely  assessments  of  the  effect  of  these  materials 
on  the  environment  and  subsequent  remedial 
cleanup  action  is  hindered  by  our  inability  to 
predict  the  forms  and  amounts  of  these  contami¬ 
nants  in  soil,  groundwater  and  plants. 

To  successfully  assess  and  repair  sites  contami¬ 
nated  with  explosives,  we  need  to  quantify  the 
mechanisms  governing  the  retention  and  trans¬ 
port  of  these  contaminants  through  the  soil  to  the 
groundwater.  Previous  investigations  of  the  fate 
of  explosive  contaminants  in  the  soil  environment 
did  not  emphasize  the  problem  of  their  potential 
mobility  to  the  groundwater.  The  ability  to  predict 
such  mobility  requires  knowledge  of  the  retention 
processes  influencing  the  fate  of  explosive  con¬ 
taminants  during  transport  in  the  soil, 

Several  studies  by  scientists  at  CRREL  have 


been  devoted  to  the  fa  te  of  explosive  contaminants 
in  soils.  These  efforts  include  the  development  of 
analytical  methodologies  tor  analyzing  explosive 
contaminants  such  as  TNT  as  well  as  DNT,  RDX 
and  HMX,  Results  of  this  work  are  available  in 
several  publications,  including  those  of  Jenkins  et 
al.  (1984,  1986,  1988)  and  Cragin  et  al.  (1985). 
Extraction  techniques  for  munitions  residues  in 
soils  have  been  recently  developed  by  Jenkins  and 
Leggett  (1985)  and  Jenkins  and  Grant  (1987), 
Palazzo  and  Leggett  (1986)  among  others  showed 
that  TNT  and  several  of  its  transformation  prod¬ 
ucts  aie  toxic  to  fish  and  several  aquatic  fauna. 

The  retention  (adsorption-desorption)  charac¬ 
teristics  of  TNT,  DNT,  RDX  and  HMX  on  pure  clay 
(Wyoming  bentonite)  were  Investigated  by  Leggett 
(1985).  He  found  that,  for  all  military  explosives 
tested,  retention  reactions  were  rapid  and  fully 
reversible.  He  examined  the  reversibility  assump¬ 
tion  by  desorption  of  the  solute  after  apparent 
equilibrium  was  attained,  In  addition,  sorption 
isotherms  appeared  to  be  linear  for  RDX  and  HMX, 
whereas  he  observed  a  moderate  nonlinearity  for 
TNT  and  DNT,  Surface  sorption  was  probably  the 
dominant  mechanism  on  the  clay  studied.  More¬ 
over,  these  contaminants  did  not  degrade  in  such 
pure  clays,  since  irreversible  behavior  was  not 
observed,  and  thus  the  total  of  the  sorbed  amounts 
was  released.  Similarly,  Penning!  on  and  Patrick 
(1990)  found  that  TNT  was  only  ,  lightly  resistant 
to  desorption,  In  fact,  after  three  sequential  des¬ 
orptions  over  12  hours,  some  88  to  93%  of  adsorbed 
TNT  was  released. 

The  above  papers  are  some  of  the  few  that 
attempted  to  quantify  the  retention  behavior  of 
TNT  and  other  explosive  contaminants  in  soils 
and  in  pure  clay .  Other  researchers  have  proposed 
mechanisms  governing  the  fate  of  TNT  and  its 
transformation  products.  Burl inson  (1 980)  showed 
that  photo-decomposition  is  the  process  forming 
1,2,5-trinitrobenzene  in  natural  waters.  Kaplan 
and  Kaplan  (1982)  proposed  microbial  transfer- 


million  as  the  dominant  mechanism  causing  the 
reduction  of  nitro  groups  to  amino  groups. 

Xhe  primary  objective  of  this  investigation  is  to 
quantify  the  potential  mobility  of  TNT  (2,4,6-trini¬ 
trotoluene)  and  RDX  (2,3,5-trinitro-l,3,5-triaz,ine) 
‘in  soils,  To  achieve  this,  we  focused  on  describing 
the  mobility  and  retention  of  TNT  and  RDX  in  a 
pure  montmorilloniteclay,  in  two  uncontaminated 
soils,  and  in  contaminated  soils  from  the  Louisiana 
Army  Ammunition  Plant,  Specifically,  we  at¬ 
tempted  to  quantify  the  retention  properties  of 
TNT  and  RDX  based  on  their  transport  (miscible 
displacement)  behavior  in  the  montmorillonite 
clay  and  soils,  We  also  tested  the  applicability  of 
mechanistic  models  that  incorporate  the  retention 
processes  with  the  convection-dispersion  trans¬ 
port  equation. 


TRANSPORT  MODELS 


equilibrium  or  a  kinetic  reaction,  where  concen¬ 
trations  in  solution  and  sorbed  phases  vary  with 
time;  reviews  of  various  forms  of  equilibrium  and 
kinetic  models  are  given  by  Murali  and  Aylmore 
(1983),  Recently,  Selim  (1992)  presented  a  review 
of  significant  features  of  sorption-exchange  reac¬ 
tions  of  the  equilibrium  and  kinetic  type,  Linear, 
Freundlich,  and  one-  and  two-site  Langmuir  equa¬ 
tions  are  most  commonly  used  to  describe equilib  • 
rium  reactions, 

The  Freundlich  equation  is  perhaps  the  sim¬ 
plest  approach  for  quantifying  the  behavior  of 
retention  of  reactive  solute  with  the  soil  matrix,  It 
is  certainly  one  of  the  oldest  of  the  nonlinear 
sorption  equations  and  has  been  used  widely  to 
describe  solute  retention  by  soils  (Helfferich  1962, 
Travis  and  Etnier  1981,  Murali  and  Aylmore  1983, 
Sposito  1984).  The  Freundlich  equation  is 

S  -  Kj  C(’  (2) 


The  classical  convective-dispersive  equation, 
which  is  generally  accepted  for  the  description  of 
dissolved  chemicals  in  the  soil  solution  (Selim 
1992),  is 


®  “  +  P  ^7  =  ®D  -  v  ~~~  -  Si  4*  i  (1) 
3<  3 1  oz 2  3z 


where 

C  =  solute  concentration  in  solution  (pg/mL) 
0  =  soil  water  content  (cm3 /cm3) 
t  -  time  (hr) 
z  =  depth  (cm) 
p  =  soil  bulk  density  (g/cm3) 

D  =  hydrodynamic  dispersion  coefficient 
(cm2 /hr) 

v  =  Darcy's  water  flux  density  (cm/hr) 

S  =  solute  concentration  associated  with  the 
solid  phase  of  the  soil  (pg/g  soil) 

4^  =  rates  of  solute  removal  (or  supply)  from  soil 
solution  (pg/cm3  •  hr)  and  not  included 
in  S, 


We  consider  the  (3S/3f)  term  as  a  fully  reversible 
process  between  the  solution  and  the  solid  phases, 
while  ‘Fj  are  rates  of  irreversible  reactions,  i.e,, 
transformations,  Processes  governing  the  interac¬ 
tions  of  individual  solute  species  need  to  be  iden¬ 
tified  if  prediction  of  the  fate  of  contaminants  in 
the  soil  using  the  convection-dispersion  equation 
(eql )  is  sought.  The  reversible  term  (3S/3f)  is  often 
used  to  describe  the  rate  of  sorption  or  exchange 
reactions  with  the  solid  matrix,  Sorption  or  ex¬ 
change  has  been  described  by  either  Instantaneous 


where 

S  =  amountofsoluteretainedby  thesoil(pg/g) 
C  =  solute  concentration  in  solution  (gg/mL) 
Kj  =  distribution  coefficient  (cm3/g) 
b  -  dimensionless  (typically  b  <  1), 

The  distribution  coefficient  describes  the  parti¬ 
tioning  of  a  solute  species  between  solid  and  liquid 
phases  over  the  concentration  range  of  interest 
and  is  analogous  to  the  equilibrium  constant  for  a 
chemical  reaction.  For  a  b  of  unity,  the  Freundlich 
equation  is  often  referred  to  as  the  linear  retention 
equation. 

Although  the  Freundlich  equation  has  been 
rigorously  derived  (Sposito  1980),  the  goodness- 
of-fit  of  the  Freundlich  equation  to  solute  retention 
data  does  not  provide  definitive  information  about 
the  actual  processes  involved,  since  the  equation  is 
capable  of  describing  data  irrespective  of  the  ac¬ 
tual  retention  mechanism.  Often  complex  reten¬ 
tion  processes  can  at  least  in  part  be  described  by 
relatively  simple  models  such  as  the  Freundlich 
equation.  Therefore,  the  Freundlich  parameters 
Kj  and  b  are  best  regarded  as  descriptive  param¬ 
eters  in  the  absence  of  independent  evidence  con¬ 
cerning  the  actual  retention  mechanism. 
Pennington  and  Patrick  (1990)  found  that  best-fit 
Freundlich  b  values  were  consistently  greater  than 
1  for  TNT  In  a  wide  range  of  soils.  Such  a  finding 
indicates  a  lack  of  sorption  maxima  for  their  range 
of  concentrations. 

The  Langmuir  isotherm  is  the  oldest  and  most 
commonly  encountered  in  soils.  It  was  developed 
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to  describe  the  adsorption  of  gases  by  solids  where 
a  finite  number  of  adsorption  sites  in  the  surface  is 
assumed  (Langmuir  1918),  As  a  result,  a  major  ad¬ 
vantage  of  the  Langmuir  equation  over  linear  and 
Freundlich  types  is  that  a  maximum  sorption  ca¬ 
pacity  is  incorporated  into  the  formulation  of 
the  model,  which  may  be  regarded  as  a  measure  of 
the  amount  of  available  retention  sites  on  the  solid 
phase.  The  standard  form  of  the  Langmuir  model  is 

=  «£..  (3) 

^max  1 

where  co  and  SmflX  are  adjustable  parameters.  Here 
to  (cm3/pg)  is  a  measure  of  the  bond  strength  of 
molecules  on  the  matrix  surface  and  Smax  (pg/g 
soil)  is  the  maximum  sorption  capacity  or  total 
amount  of  available  sites  per  unit  soil  mass,  In  an 
attempt  to  classify  the  various  shapes  of  sorption 
isotherms,  Sposito  (1984)  recognized  that  the 
Langmuir  isotherm,  or  L-curve  isotherm,  is  the 
most  commonly  used.  The  Langmuir  sorption  iso¬ 
therm  has  been  used  extensively  by  scientists  for 
several  decades.  Travis  and  Etnier  (1981)  pro¬ 
vided  a  review  of  studies  where  the  Langmuir 
isotherm  was  used  to  describe  P  retention  for  a 
wide  range  of  soils.  Moreover,  Langmuir  isotherms 
were  used  successfully  to  describe  Cd,  Cu,  Pb  and 
Zn  retention  in  soils,  Leggett  (1985)  and  Pennington 
and  Patrick  (1990)  successfully  used  the  Langmuir 
approach  to  describe  batch  results  of  TNT  reten¬ 
tion  on  a  reference  clay  and  on  soils  having  a  wide 
range  of  properties. 

Multisite  or  multireaction  models  deal  with  the 
multiple  interactions  of  one  species  in  the  soil 
environment.  Such  models  are  empirical  and  are 
based  on  the  assumption  that  a  fraction  of  the  total 
sites  is  highly  kinetic,  whereas  the  remaining  sites 
interact  slowly  or  instantaneously  with  the  chemi¬ 
cals  in  the  soil  solution  (Selim  et  al.  1976,  Jardine  et 
al.  1985).  Nonlinear  equilibrium  (Freundlich)  was 
assumed  to  describe  retention  reactions  associ¬ 
ated  with  equilibrium  sites,  whereas  first-  or  nth- 
order  kinetic  reactions  were  the  mechanisms  associ¬ 
ated  with  the  kinetic  sites,  Such  a  two-site  ap¬ 
proach  proved  successful  In  describing  observed 
extensive  tailing  of  breakthrough  results. 

Another  two-site  approach  was  proposed  by 
Theiset  al.  (1988)  for  Cd  mobility  and  adsorption 
on  geothite.  They  assumed  that  the  nature  of  reac¬ 
tions  is  governed  by  second-order  kinetic  reactions, 
The  reactions  were  assumed  to  be  consecutive 
where  the  second  reaction  was  irreversible, 
Amacher  et  al.  (1988)  developed  a  multireaction 


model  that  includes  concurrent  and  concurrent- 
consecutive  processes  of  the  noniinoa  r  kinetic  type. 
The  modei  was  capable  of  describing  the  retention 
behavior  of  Cd  and  Cr(VI)  with  time  for  several 
soils.  In  addition,  the  model  predicted  that  a  frac¬ 
tion  of  these  heavy  metals  was  irreversibly  retained 
by  the  soil.  Recently,  Amacher  et  al.  (1990)  con¬ 
cluded  that  the  multireaction  model  could  also  suc¬ 
cessfully  describe  adsorption  of  Hg  for  several  soils. 

In  the  Multi-Reaction  and  Transport  Model 
(MRTM)  of  Selim  et  al.  (1990),  a  solute  species  was 
assumed  to  be  presen  t  in  the  soil  solution  phase  (C) 
and  in  four  phases  representing  solute  retained  by 
the  soil  matrix  ns  Se,  Sj,  S2  and  S'lrr,  The  model 
further  assumed  that  Se,  Sj  and  S2  were  in  direct 
contact  with  the  solution  phase  and  were  gov¬ 
erned  by  concurrent  type  reactions.  Here,  Se  repre¬ 
sents  the  amount  of  solute  that  is  sorbed  reversibly 
and  is  in  equilibrium  with  C  at  all  times.  The 
governing  equilibrium  retention-release  mecha¬ 
nism  was  that  of  the  nonlinear  Freundlich  type  as 
discussed  previously. 

The  retention-release  reactions  associated  with 
S-,  and  S2  had  to  be  In  direct  contact  with  C  to  be 
governed  by  reversible  processes  of  the  (nonlin- 


ear)  kinetic  type 

Se  =  Kd  Cb 

(4) 

^-L  shSc*  -  fc2  Si 

(5) 

at  p 

^2=  h@-cm  -fc4  s2 
at  '  p 

(6) 

where  fcj  to  k4  are  the  associated  rate  coefficients 
(hr-1).  These  two  phases  (S^  and  S2)  may  be  re¬ 
garded  as  the  amounts  sorbed  on  surfaces  of  soil 
particles  and  chemically  bound  to  Al  and  Fe  oxide 
surfaces  or  other  types  of  surfaces,  although  it  is 
not  necessary  to  have  a  priori  knowledge  of  the 
exact  retention  mechanisms  for  these  reactions  to 
be  applicable.  Moreover,  these  phases  may  be 
characterized  by  their  sorption  and  release  kinet¬ 
ics  to  the  soil  solution  and  thus  are  susceptible  to 
leaching  in  the  soil.  In  addition,  primary  differ¬ 
ences  between  these  two  phases  not  only  lie  in  the 
difference  in  their  kinetic  behavior  but  also  in  the 
degree  of  nonlinearity,  as  indicated  by  the  param¬ 
eters  it  and  in,  The  multireaction  model  also  con¬ 
siders  irreversible  solute  removal  via  a  retention 
sink  term  Q  to  account  for  irreversible  reactions 
such  as  precipitation-dissolution  and  mineraliza¬ 
tion-immobilization,  among  others.  Weexprossed 
the  sink  term  as  a  first-order  kinetic  process 


3 


Q  =  P  ^  =  *»  ©c  (7) 

where  is  the  associated  rate  coefficient  (hr1). 

We  incorporated  the  above  retention  mecha¬ 
nisms  into  the  classical  convection-dispersion 
equation  (eq  1)  to  predict  solute  retention  as  gov¬ 
erned  by  the  multireaction  model  during  trans¬ 
port  in  soils  (Selim  et  al,  1989).  The  initial  and 
boundary  conditions  used  were 

C  =  0,  f  =  0,  OsSzsiL  (8) 

se  =  Si  =  S2  =  o,  t  =  0,  0  ^  z  *  L  (9) 


oC0  =  -D  —  +  vC,  z  =  0,  t  <  T  (10) 

dz 

0  =  -D^-  +  vC,  z  =  0,  fisT  (11) 

dz 

3C.0,  z  -  L,  tao  (12) 

dz 

where  it  is  assumed  that  a  solute  solution  of  a 
concentration  (C0)  is  applied  to  a  soil  column 
having  a  length  L  for  a  given  duration  T  and  was 
thereafter  followed  by  a  solute-free  solution.  The 
numerical  sol  utlon  for  the  multireaction  and  trans¬ 
port  model  is  documented  in  the  program  MRTM 
and  is  given  in  Selim  et  al.  (1990). 


EXPERIMENTAL 
Soils  and  reagents 

A  contaminated  soil  was  collected  from  the 
Louisiana  Army  Ammunition  Plant  in  north  Loui¬ 
siana  .  It  appeared  to  be  disturbed  and  m  ixed;  thus, 
it  was  not  possible  to  exactly  identify  soil  layers. 
An  uncontaminated  (original)  soil  was  also 
sampled  about  300  m  from  the  contaminated  site, 
The  surface  (0-30  cm)  horizon  of  the  uncontami¬ 
nated  (original)  soil  was  sampled  to  serve  as  a 
blank  (a  check)  for  soil  physical  and  chemical 
properties  in  comparison  to  the  contaminated 
sample.  Two  subsamples  from  the  contaminated 
and  uncontaminated  soils  were  forwarded  to 
CRREL  for  complete  analysis  of  total  TNT  and 
RDX  contents  as  well  as  other  explosive  contami¬ 
nants.  The  contaminant  analyses  revealed  that  no 
TNT,  RDX  or  other  contaminants  were  present  in 
the  uncontaminated  soil,  The  soil  at  the  site  is 
Kolin  soil  (fine-silty,  siliceous,  Thermic  Glossaqulc 
Paleudalfs),  Further  details  may  be  obtained  from 


the  document  entitled  Soil  Survey  of  Caddo  Parish, 
Louisiana  (USDA-SCS 1980)  and  from  Amacher  et 
al.  (1989). 

A  reference  clay  was  also  used  in  this  study: 
SWy-lmontmorillonite(CrookCounty,  Wyoming, 
bentonite)  obtained  from  theSource  Clays  Reposi¬ 
tory,  University  of  Missouri,  Sposito  et  al.  (1983a,  b) 
and  Gaston  and  Selim  (1990)  have  presented  ex¬ 
tensive  data  for  ion  exchange  (of  Ca,  Mg  and  Na) 
on  pure  montmorillonite  fractionated  from  this 
material,  Most  significant  is  Leggett's  (1985)  ex¬ 
tensive  work  on  theadsorption-desorption  of  TNT 
and  RDX  on  pure  montmorillonite  clay,  which 
indicated  fully  reversible  TNT  and  RDX  retention 
mechanisms,  Based  on  this  knowledge  ofTNT and 
RDX  in  this  reference  clay,  we  conducted  trans¬ 
port  and  adsorption-desorption  experiments  to 
determine  the  mobility  of  these  contaminants  in 
soils.  Our  primary  objective  was  to  identify  the 
sorption-release  mechanisms  of  TNT  on  2.1  mont- 
morillonitic  clays.  To  predict  the  potential  mobil¬ 
ity  of  TNT  and  RDX  in  soils  under  field  conditions, 
we  also  evaluated  the  extent  of  time-dependent  or 
kinetic  behavior  of  the  retention  mechanisms. 
Mineral  compositions  of  the  less  than  2-pm  frac¬ 
tions  of  the  reference  clay  material  were  deter¬ 
mined  from  X-ray  diffractograms  of  oriented  slide 
mounts  after  K-saturation  with  heat  treatments  at 
300  and  550°C  and  Mg-saturation  with  glycolation, 
X-ray  analyses  were  run  from  2-32°  20,  using  Cu- 
Ka  radiation  with  a  vertical  goniometer  equipped 
with  a  20  compensator  slit  and  graphite 
monochronometer,  Estimates  of  composition  were 
based  on  areasof  characteristic  diffraction  peaks — 
it  was  99%  montmorillonite. 

A  Norwood  soil  was  also  chosen  for  this  study 
for  several  reasons.  This  soil  is  predominant  in  the 
Red  River  valley,  which  is  adjacent  to  the  contami¬ 
nated  site  near  the  Louisiana  Army  Ammunition 
Plant,  Norwood  soil  isclassified  as  fine  silty,  mixed, 
calcareous,  Thermic  and  Typie  Udifluvent  with 
79%  sand,  18%  silt  and  3%  clay .  In  addition,  previ¬ 
ous  work  on  this  soil  indicates  that  it  is  of  low 
organic  matter  content  (0,32%),  high  pH  (7.4)  and 
low  cation  exchange  capacity  (4.1  cmolc/kg),  and 
it  is  highly  permeable  to  water  infiltration  with  a 
low  water  retention  capacity  (Amacher  et  al.  1988, 
Buchter  et  al,  1989).  The  adsorption  behaviors  of 
several  heavy  metals  in  this  soil  have  been  recently 
documented  by  Buchter  et  al.  (1989),  This  soil  is  a 
worst  case,  owing  to  its  high  leaching  capability 
and  thus  the  potential  for  contamination  of  the 
groundwater, 

Analytical  standards  for  TNT  were  prepared 
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from  a  reagent  chemical  purchased  from  Eastman 
Kodak.  The  RDX  used  was  provided  by  CRREL 
scientists.  Both  RDX  and  TNT  were  dried  to  a 
constant  weight  over  dry  calcium  chloride  in  a 
desiccator  in  the  dark.  Standards  were  prepared  in 
methanol  for  TNT  and  acetonitrile  and  methanol 
for  RDX.  A  combined  analyte  calibration  standard 
was  prepared  in  methanol.  Volumetric  flasks  were 
stoppered,  and  the  flask  closures  were  wrapped 
with  parafilm  to  minimize  evaporation  of  solvent. 
All  standards  were  stored  in  the  dark  at  4°C. 

All  analyses  for  TNT  and  RDX  in  soil  water 
samples  were  carried  out  using  HPLC  methodol¬ 
ogy  as  established  by  CRREL.  Briefly,  all  HPLC 
analyses  were  conducted  using  an  IBM  LC/9533 
Ternary  Gradient  Liquid  Chromograph  with  an 
IBM  LC/9523  variable-wavelength  UV  detector 
set  at  a  wavelength  (X)  of 244  nm.  A  Rheodyne  7125 
injector  valve  with  a  100-p.Lsample  loop  was  used 
and  samples  were  introduced  to  the  column  by 
overfilling  the  loop.  A  Supelco  25-cm  x  4.6-mm  i.d. 
(5-pm)  reversed-  phase  column  (LC-18)  was  used 
with  an  eluent  of  methanol  and  water  (50:50)  at  a 
flow  rate  of  1.5  mL/min, 

A  series  of  standards  was  prepared  in  methanol 
from  the  combined  analyte  calibration  standard. 
The  concentrations  included  0.5, 1, 2, 5, 10  and  20 
mg/L  of  RDX  and  TNT.  Initial  calibration  was 
done  by  analyzing  these  concentrations  in  du¬ 
plicate  and  in  random  order.  The  calibration  curve 
was  determined  to  be  linear,  with  the  slope  of  the 
regression  line  the  response  factor.  Final  calibra¬ 
tion  was  done  using  a  single-point  working  stan¬ 
dard,  which  was  analyzed  in  replicate  at  the  begin¬ 
ning,  the  middle  and  the  end  of  analysis  daily.  The 
mean  peak  height  was  used  to  calculate  the  re¬ 
sponse  factor,  which  was  compared  to  the  re¬ 
sponse  factor  from  the  initial  calibration. 

Sorption  experiments 

The  retention  of  applied  TNT  during  sorption 
was  investigated  using  the  batch  method  as  out¬ 
lined  by  Amacher  et  al.  (1988),  Initial  concentra¬ 
tions  of  TNT  in  the  solution  added  to  the  soil  were 
1,110,  2.275,  5.559,  7,708,  10.22  and  11,31  mg/L. 
Background  solution  composition  was  0.005-M 
Ca(NC>3)2.  For  the  SWy-1  montmorillonite,  TNT 
retention  was  carried  out  as  folio  ws.  Triplicate3.5- 
g  samples,  consisting  of  0,5-g  of  SWy-1  nvmtmo- 
rillonite  and  3  0  g  of  acid-washed  sand,  were  added 
to  40-mL  polypropylene  centrifuge  tubes  of  known 
weight,  The  acid-washed  sand  was  assumed  to  be 
nonreactive,  which  was  necessary  to  maintain 
mixing  and  distribution  of  the  SWy-1  rnontmoril- 


lonite  (see  Gaston  and  Selim  1990).  Twenty-five 
milliliters  of  TNT  solution  were  added  to  the  trip¬ 
licate  samples,  which  were  vortex  mixed.  The 
tubes  were  placed  endwise  on  a  shaker  set  to  sha  ke 
at  120  oscillations/minute.  The  samples  were 
shaken  for  15  minutes  every  2  hours.  After  3  days, 
the  samples  were  centrifuged  and  the  supernatant 
decanted.  Concentrations  of  TNT  in  the  sample 
solution  were  based  on  the  HPLC  method  de¬ 
scribed  above.  Amounts  of  the  TNT  sorbed  or 
retained  by  the  SWy-1  montmorillonite  samples 
were  calculated  from,  the  difference  between  ini¬ 
tial  (applied)  concentrations  of  each  sample, 
Correction  was  also  made  for  the  entrained  super¬ 
natant. 

Desorption  experiments 

Desorption  or  release  of  TNT  was  carried  out 
following  sorption  as  outlined  by  Amacher  et  al, 
(1988),  with  successive  dilution  steps,  After  the 
solution  was  decanted,  the  first  dilution  step  was 
carried  out  by  adding  25  mL  of  Ca(N03)2  to  each 
test  tube,  The  tubes  were  then  vortex  mixed  and 
placed  on  the  shaker  as  described  in  the  above 
section,  The  samples  were  shaken  for  15  minutes 
every  2  hours.  After  a  minimum  release  time  of  4 
to  6  hours  for  each  step,  the  samples  were  centri¬ 
fuged  and  the  supernatant  decanted  for  TNT  analy¬ 
sis.  Five  dilution  steps  were  carried  out  using  this 
procedure.  In  addition  to  the  TNT  concentrations 
after  each  release  or  extraction  step,  the  amount  of 
TNT  released  and  that  remaining  on  matrix  sur¬ 
faces  were  calculated  on  the  basis  of  solution  con¬ 
centrations  and  input  (initial)  amounts  added. 
Release  beyond  the  fifth  dilution  step  was  not 
done  because  the  TNT  in  the  soil  solution  dropped 
below  detection  limits. 

For  the  contaminated  soil  (samples  A  and  B), 
the  above  method  was  followed  except  that  3-g 
soil  samples  were  used.  Again,  triplicate  soil  sam¬ 
ples  were  used  and  the  desorption  reaction  dura¬ 
tion  was  1  day  for  each  extraction  step,  after  which 
solutions  were  analyzed  for  TNT  and  RDX.  De¬ 
sorption  was  carried  out  for  8  days  for  eight  extrac¬ 
tions.  Release  reaction  was  terminated  after  the 
eighth  day,  as  TNT  and  RDX  levels  in  the  extrac¬ 
tion  solutions  were  below  the  HPLC  detection 
limit. 

Transport  experiments 

Transport  experiments  were  conducted  to  de¬ 
scribe  the  potential  mobility  of 'ENT  and  RDX  in  a 
standard  reference  material  (pure  clay),  in  con¬ 
taminated  and  uncontaminated  soils  obtained  from 
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Table  1.  Parameters  of  the  various  soil  columns  for  the  TNT  and  RDX  miscible 
displacement  experiments, 


Soil 

Experiment 

number 

WtUer  content  0 
(ew-Vcm-V 

Hulk  density  p 
(g/em3) 

Flux  v 
(em/hr) 

Bentonite-sand 

102 

0.363 

1 ,56 

0.294 

Bentonite-sand 

110 

0.385 

1.56 

0.295 

Bentonite-sand 

111 

0.385 

1.56 

0.295 

Norwood 

103 

0.403 

1.46 

0.900 

Norwood 

105 

0.404 

1.47 

0,747 

Kalin* 

112 

0.466 

1.45 

0.302 

Kolln 

sample  A** 

108 

0.510 

1,42 

0,279 

Kolln 

sample  B*‘ 

107 

0.517 

1.45 

0.310 

*  Uncontaminated  sample  from  Louisiana  Army  ammunition  site. 
**  Contaminated  samples  from  Louisiana  Army  ammunition  site. 


Table  2.  Concentration  of  TNT  and  RDX  input  pulses,  background  solution  and  pulse 
duration  for  the  miscible  displacement  column  experiments. 


Soli 

Experiment 

number 

Input 

pulse 

Concent¬ 

ration 

(Mg/mU 

Background 

solution 

Input 

pulse 

(pore  volumes) 

Duration 

(hr) 

Bentonite- 

sand 

102 

TNT 

100.00 

50% 

methanol 

4.33 

54.6 

Bentonite- 

sand 

110 

TNT 

10,65 

0.005-M 

Ca(NO,)2 

7.51 

95.9 

Bentonite- 

sand 

111 

RDX 

10.65 

0.005-M 

Ca(N03)2 

4.33 

90.1 

Norwood 

103 

TNT 

100.00 

0.005-M 

Ca(NO,i)j 

7,11 

30.6 

Norwood 

105 

TNT 

10,28 

0.005-M 

Ca(N03)2 

5.18 

28.1 

Kolln* 

112 

TNT 

10.35 

0,005-M 

Ca(NOj)j 

5,00 

74.5 

RDX 

10.04 

0.005-M 

Ca(NOj)j 

5.00 

74.5 

Kolln 
sample  A 

108 

— 

— 

0,005-M 

Ca(NOj)j 

37.9 

693 

Kolln 
sample  B  * 

107 

-- 

— 

0.005-M 

Ca(NO|)a 

40.5 

678 

‘  Uncontammatcd  sample  from  Louisiana  Army  Ammunition  Plant  site.  Applied  pulse  contained  both  RDX  and  TNT, 
**  Contaminated  samples  from  Louisiana  Army  Ammunition  Plant  site. 
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the  Louisiana  Army  Ammunition  Plant,  and  in  the 
highly  pormoablo  Norwood  soil.  For  the  standard 
clay,  wo  used  a  1:9  bonlonito-sand  mixture.  A 
mixture  with  acid-washed  sand  assured  unifor¬ 
mity  and  allowed  water  flow  through  the  packed 
soil  columns,  Our  transport  experiments  used  the 
miscible  displacement  approach  described  by 
Selim  et  al,  (1989),  where  leaching  through  packed 
soil  columns  was  monitored.  Each  soil  was  air- 
dried  and  passed  through  a  2-mm  sieve  before  use. 
Plexiglas  columns  (6,4  cm  in  diameter  x  10  cm  in 
length)  were  uniformly  packed  witheachsoil  to  an 
appropriate  bulk  density .  Constant  flux  was  main¬ 
tained  with  variable-speed  piston  pumps  (fluid 
metering  model  RPG6),  A  fraction  collector  was 
used  to  collect  column  effluent  (ISCO  model  68- 
2180)  at  different  intervals  that  depended  on  efflu¬ 
ent  fluxes. 

Initially,  soil  columns  were  slowly  wet  with 
0,005-M  Ca(N03)2  and  leached  with  this  solution 
until  negligible  Mg,  K  or  Na  appeared  in  the 
effluent.  Analysis  of  leachate  collected  in  the  Ca- 
saturation  step  indicated  that  essentially  all  ini¬ 
tially  extractable  Mg,  K  and  Na  had  been  removed, 
Over  20  pore  volumes  of  0,005-M  Ca(N03)2  were 
introduced  into  each  column  at  a  constant  flux  to 
equilibrate  the  soils  with  the  background  solution. 
Comparison  ol  column  weights  after  wetting  with 
Ca  solutions  to  initial  air-dry  weights  indicated 
that  nearly  full  saturation  had  been  achieved.  Ex¬ 
perimental  parameters  of  the  column  leaching 
studies  are  given  in  Tables  1  and  2. 

A  pulse  of TNT  having  a  concentration  of  either 
10  or  100  mg/L  in  a  0.005-M  Ca(N03)2  solution 
was  introduced  into  each  column  at  a  constant 
flux.  The  applied  TNT  was  followed  by  several 
pore  volumes  of  the  background  solution  of  0.005- 
M  Ca(N03)2  to  elute  the  columns  and  obtain  a 
complete  TNT  breakthrough  curve  for  each  soil. 
The  experimental  miscible  displacement  proce¬ 
dure  just  described  was  followed  for  all  packed 
soil  columns  and  the  pure  clay  with  only  one 
exception:  forsoil  columns  with  thecontaminated 
samples  A  and  13,  pulse  applications  of  TNT  were 
not  made.  Instead,  we  monitored  the  leaching  of 
TNT  and  RDX  that  was  already  present  in  the 
contaminated  soils.  In  addition,  the  backgrourd 
solution  of  0.005-M  Ca(N03)2  was  introduced  to 
the  air-dry  soil  and  the  extent  of  leaching  or  release 
based  on  concentrations  of  TNT  and  RDX  in  the 
effluent  solution  was  monitored. 

To  determine  the  hydrodynamic  dispersion 
coefficient  (D)  of  eq  1,  a  pulse  of  tritium  (-1H20) 
was  applicu  to  the  reference  clay  columns  prior  to 


the  TNT  pulses.  The  pulse  volume  was  equivalent 
to  one  pore  volume  and  the  background  solution 
was  0.005-M  Ca(NQ3)2,  which  was  maintained  at 
all  times.  Liquid  scintillation  spectrometry 
(Beckman  LS7000)  was  used  to  count  the  tritium 
radio  tracers.  Best-fit  model  parameters  for  D  were 
obtained  by  use  of  nonlinear  least-squares  optimi¬ 
zation  methods  of  the  tritium  results  (see  van 
(Jenuchten  1981).  The  estimate  for  D  was  0,00134 
m2  day-1,  and  that  for  a  retardation  factor  was 
close  to  unity. 

RESULTS  AND  DISCUSSION 
Sorption  isotherms 

Figure  1  is  the  sorption  isotherm  for  TNT  on  the 
SWy-1  montmorillonite.  This  figure  represents  the 
amount  of  TNT  sorbed  vs.  that  in  the  soil  solution 
following  1  day  of  retention  reaction.  The  amount 
retained  by  the  clay  material  is  expressed  in  terms 
of  micrograms  of  TNT  retained  per  gram  of  SWy- 
1  montmorillonite  (|ig/g).  The  acid-washed  sand 
was  considered  nonreactive  and  was  thus  not 
included  in  any  of  the  calculations  reported  below. 

Isotherm  results  of  Figure  1  were  described  in 
three  different  ways.  First  a  linear  Freundlich  form 
was  used 

S  =  KdC  (13) 

where 

S  =  amount  sorbed  (pg/g) 

C  =  solution  concentration  (pg/mL) 

Kd  =  distribution  coefficient  (em3/g), 

The  linear  form  is  similar  to  the  Freundlich  equa¬ 
tion  (eq  2)  with  b  =  1.  The  Kj  parameter  describes 
the  partitioning  of  TNT  between  the  solid  and 
liquid  phases  over  the  concentration  range  and  is 
analogous  to  the  equilibrium  constant  for  a  chemi¬ 
cal  reaction,  The  best  fit  Kd,  which  provided  the 
solid  line  in  Figure  1,  was  51.722  ±  2.929  cnvVg, 
There  is  strong  deviation  from  linear  behavior, 
especially  for  low  TNT  concentrations  (less  than  5 
pg/mL)  as  depicted  in  Figure  1.  As  a  result,  the 
common  nonlinear  Freundlich  form  was  used 
(Selim  et  al.  1990),  Nonlinear  least-square  optimi¬ 
zation  was  carried  out  to  provide  estimates  for  Kj 
and  b  based  on  eq  2.  Parameter  values  that  pro¬ 
vided  the  best  fit  for  the  isotherm  results  were 
of  65.881  ±  11.046  cmVg  and  h  of  0.841  ±0,101, 
Isotherm  prediction  using  the  nonlinear 
Freundlich  form  is  shown  by  the  solid  cun  e  of 
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TNT  Solution  Concentration  (ng/mL) 


Figure  1.  TNT  sorption  on  bentonite  clay  (SWy-1  montmorillonite,  Crook  County ,  Wyo¬ 
ming).  The  solid  line  is  a  fitted  linear  model. 


TNT  Solution  Concentration  (ng/mL) 


Figure  2  TNT  sorption  on  bentonite  clay  (closed  circles)  and  fitted  Freundlich  and  Langmuir 
models. 
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Figure  3,  TNT  sorption  an  bentonite  clay  (closed  circles)  and  fitted  modified  Langmuir  model. 


Figure  2.  Such  a  nonlinear  approach  provided 
good  predictions  for  the  low  TNT  concentration 
range .  Greater  discrepancy  between  predicted  and 
experimental  results  were  obtained  as  TNT  con¬ 
centrations  increased.  The  b  value  of  0.84  for  the 
nonlinear  form  is  consistent  with  isotherms  for 
other  dissolved  chemicals  for  other  soils  and  pure 
minerals.  In  general,  the  parameter  b  is  typically 
less  than  unity,  as  was  found  in  this  study  (Buchter 
et  al,  1989).  In  contrast,  Pennington  and  Patrick 
(1990)  found  that  best-fit  Freundlich  b  values  for 
TNT  sorption  were  consistently  greater  than  unity 
for  a  wide  range  of  soils. 

Langmuir  predictions  were  also  used  to  de¬ 
scribe  the  isotherm  data  for  TNT  sorption  on  the 
clay.  The  best-fit  prediction  is  indicated  by  the 
dashed  curve  in  Figure  2,  It  is  obvious  that  this 
approach  is  inappropriate  for  describing  the  iso¬ 
therm  results  in  the  concentration  range  shown, 
The  values  for  co  and  Smax,  which  provided  the 
best  fit  shown,  are  0,196  ±  0.036  (mL/pg)  and 
493.467  ±47.905  (pg/g)  respectively.  Leggett  (1985) 
found  that  the  use  of  an  extended  Langmuir  form 
as 


5  _  (oC  +  k  C 

^max  (1+W  C) 


(14) 


provided  the  best  predictions  of  TNT,  DNT,  RDX 
and  HMX  isotherms  for  a  similar  concentration 
range  (tc  is  an  arbitrary  parameter).  The  above 
formulation  may  be  regarded  as  a  two-site  ap¬ 
proach,  where  a  combination  of  linear  plus 
Langmuir  processes  describes  the  retention  mecha¬ 
nisms.  For  our  TNT  sorption  isotherm,  the  use  of 
the  combined  approach  improved  the  prediction 
capability  of  the  experimental  data,  as  shown  in 
Figure  3.  Parameter  estimates  that  provided  best- 
fit  predictions  of  Figure  3  were  37.258  ±  22.698, 
1.502  ±  1.189  and  44.477  ±  4.681  for  Smax,  m  and  k 
respectively .  In  fact,  the  predicted  curve  was  closely 
similar  to  that  for  the  Freundlich  approach.  How¬ 
ever,  the  large  values  for  the  parameter  standard 
errors  given  in  Table  3  arc  indications  that  the 
modified  Langmuir  approach  provides  inadequate 
representation  of  the  TNT  retention  mechanisms 
on  the  reference  clay. 

Desorption  and  hysteresis 

Following  the  sorption  step,  TNT  was  des¬ 
orbed  using  ;. .?  successive  dilution  method  as 
outlined  by  Amacher  et  al.  (1988).  Desorption 
isotherms  for  two  initial  applied  concentrations 
are  shown  in  Figure  4.  Each  desorption  curve  is  a 
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Table  3.  Best-fit  model  parameters  for  batch  sorption  data  of  TNT  on 
reference  clay  (SWy-1  montmorillonite). 


Model 

Rate 

coefficient 

Parameter 

estimate 

error 

r^ 

Linear 

Kd  (cm3/ g) 

51.722 

1.197 

0.995 

Freundlich 

Kd  (cmtyg) 

65,881 

4,297 

0,996 

b 

0,841 

0.039 

Langmuir 

^max 

493  46 

19.579 

0.976 

(0  (mL/gg) 

0.917 

0.0146 

Modified 

(MR/?) 

37.258 

8,8307 

0.998 

Langmuir 

to  (mL/ng) 

1.502 

0,4624 

k  (mL/gg) 

44.477 

1,8209 

result  of  several  dilution  steps,  each  followed  by 
equilibration  for  4  to  6  hours.  Another  way  of 
representing  the  desorption  results  is  in  terms  of 
concentration  of  TNT  in  the  soil  solution  and  in  the 
sorbed  phase  vs.  the  number  of  extractions  (i.e,, 
the  desorption  step),  These  results  are  shown  in 
Figures  5  and  6  respectively.  As  extraction  contin¬ 
ued,  the  amount  released  for  the  various  TNT 
initial  concentrations — i.e.,  the  amount  desorbed 
to  the  soil  solution— continued  to  decrease,  The 


initial  TNT  concentrations  ranged  from  1.11  to 
11.31  |ig/mL.  In  the  meantime,  owing  to  the  suc¬ 
cessive  dilutions,  TNT  concentrations  and  the  to¬ 
tal  TNT  content  retained  by  the  clay  continued  to 
decrease. 

A  mass  balance  was  attempted  for  each  batch 
experiment  and  for  all  TNT  initial  concentrations 
using  the  total  amount  released  vs.  that  applied. 
These  mass-balance  calculations  indicated  that 
TNT  retention  is  fully  reversible,  Thus,  for  the 


10 


(fig/g  soil)  <«  ■*,  TNT  in  Solution  (jig/mL) 


Number  of  Extractions 


gure  5.  TNT  concentration  in  solution  in  the  bentonite  clay  vs,  number  of  dilutions- 


iruLiiuriz  uunrig  Mcai/rpwvrii 


Number  of  Extractions 

Figurcb,  Amount  of  TNT  sorbed  on  bentonite  clay  vs,  number  of 'dilutions-extractions  during 
desorption. 
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TNT  Solution  Concentration  (pg/mL) 


Figure  7.  TNT  sorption-desorption  isotherms  for  bentonite  clay  for  all  initial  TNT  con¬ 
centrations  used, 


bentonite  clay  used  here,  there  were  no  apparent 
irreversible  interactions  of  TNI  each  as  degrada¬ 
tions.  However,  this  is  not  the  case  for 
inhomogeneous  soils,  where  organic  matter  at 
various  levels  of  decomposition  and  various  clay 
minerals  and  oxides  play  significant  roles  in  the 
chemical,  physical  and  biological  Interactions 
within  the  soil  environment.  The  level  of  biologi¬ 
cal  activities  in  nonhomogeneous  soils  is  demon¬ 
strated  by  the  extent  of  degradation  potential  of 
military  explosives  (fora  review  see  final  report  by 
Weston  1988). 

Since  retained  TNT  was  completely  desorbed 
during  a  reaction  time  of  only  2  days,  one  may 
regard  the  retention  process  as  fast  or  instanta¬ 
neous,  where  local  equilibrium  is  reached  in  a 
relatively  short  reaction  time  (Rubin  1983),  To 
further  test  the  hypothesis  that  the  Local  Equilib¬ 
rium  Assumption  (LEA)  is  valid  for  TNT  retention 
with  the  bentonite  clay,  one  may  examine  the 
sorption  and  desorption  isotherms  shown  in  Fig¬ 
ure  4.  The  fact  that  the  desorption  isotherms  do  not 
deviate  appreciably  from  the  sorption  isotherm  is 
a  definite  indication  of  reversibility  as  well  aB  of 
local  equilibrium  conditions  as  the  dominant  re¬ 
tention  mechanism,  The  nonsingularity  orhyster- 
etic  behavior  of  the  TNT  Isotherms  may  be  a  result 


of  kinetic  behavior  of  reversible  reaction  mecha¬ 
nisms,  irreversible  mechanisms  such  as  degrada¬ 
tion,  or  artifacts  from  the  technique  used  for  the 
desorption-release  process.  Selim  et  al.  (1976)  have 
shown  the  effect  of  kinetic  behavior  on  the  ob- 
served  nonsingularity  of  sorption-desorption  iso¬ 
therms,  The  TNT  results  for  all  desorption  iso¬ 
therms  are  given  in  Figure  7  and  indicate  limited 
deviations  between  the  sorption  and  desorption 
curves  for  all  initial  TNT  concentrations  investi¬ 
gated. 

Transport 

The  transport  of  TNT  in  miscible  displacement 
columns  packed  with  a  bentonite-sand  mixture  is 
illustrated  by  the  results  shown  in  Figure  8.  In  this 
case,  we  introduced  into  the  bentonite-sand  col¬ 
umn  a  pulse  of  TNT  having  a  concentration  (C0)  of 
100  mg/Lin  a  background  solution  of  50%  metha¬ 
nol  in  distilled  water,  The  volume  of  the  pulse 
input  and  other  characteristics  of  the  miscible 
displacement  column  are  given  in  Table  2,  In  Fig¬ 
ure  8,  the  results  are  presented  in  terms  of  the 
relative  TNT  concentration  (C  /  C0)  versus  number 
of  pore  volumes  ( V /  V0),  where  is  the  volume  of 
pore  space  in  the  column, 

The  results  indicate  the  early  arrival  of  TNT  in 


12 


Pore  Volume  (V/V0) 

Figure  8.  TNT  breakthrough  results  (dosed  circles)  fora  bentonite-sand  column  with  input  TNT 
concentration  (CJ  of  WO  pg/mL  in  a  50%  methanol  background  solution.  The  solid  curve  is  a 
model  prediction  using  the  transport  equation  with  a  linear  retention  model. 


the  effluent  solution.  Moreover,  the  peak  concen¬ 
tration  exceeded  95%  of  the  input  concentration 
(C0).  The  TNTmobility  for  this  case  resembles  that 
of  a  nonreaclive  solute,  where  the  recovery  of 
applied  TNT  in  the  effluent  solution  exceeded 
95%.  The  amount  of  TNT  recovery  was  estimated 
based  on  the  area  under  the  breakthrough  results, 
Very  little  retardation  of  the  TNT  pulse  was  ob¬ 
served  here.  This  is  confirmed  by  the  relative  TNT 
concentration  in  the  effluent  (C/C0)  reaching  0.50 
after  only  1,7  pore  volumes  (V/V0).  For  a 
nonreactlve  solute,  such  relative  concentration  is 
achieved  at  approximately  V/V0  of  unity.  Such 
results  directly  indicate  the  effect  of  methanol  as  a 
solvent  on  the  extent  of  reaction  of  TNT  in  the 
bentonite-sand  column. 

Attempting  to  describe  the  transport  behavior 
of  TNT  in  the  bentonite-sand  column  when  metha¬ 
nol  was  present  (Fig.  8),  we  found  that  a  simple 
linear  retention  model  provided  a  good  descrip¬ 
tion  of  the  breakthrough  results.  The  K j  value  that 
gave  the  best  fit  to  the  data  was  0.2036  ±  0,0193 
cm3/g,  which  indicates  a  low  parameter  standard 
of  error  for  the  estimated  Kj.  The  goodies  of  fit  of 
the  linear  model  to  the  experimental  data  was 
based  on  r2  values  as  calculated  us'ng  the  principle 


of  extra  sum  of  squares  (Amacher  et  al.  1988),  The 
r2  value  for  the  results  of  Figure  8  was  0.987.  We 
should  point  out  that  the  linear  model  could  not 
describe  the  slow  release  of  TNT,  as  indicated  by 
the  tailing  of  the  breakthrough  results.  It  is  con¬ 
ceivable  that  a  small  proportion  of  the  desorbed 
TNT  was  strongly  sorbed  by  the  clay  matrix  sur¬ 
faces,  Such  strongly  held  TNT  forms  may  account 
for  the  slow  release  during  leaching  in  the  bento¬ 
nite-sand  column, 

The  TNT  breakthrough  results  shown  in  Figure 
9a  reflect  the  mobility  of  TNT  where  the  applied 
pulse  concentration  (C0)  was  10.65  mg/L  and  the 
background  solution  was0.005-MCa(N03)2.  Fig¬ 
ure  9a  shows  that  the  TNT  movement  was  re¬ 
tarded:  arrival  in  the  exit  solution  was  not  detected 
until  after  seven  pore  volumes,  and  peak  concen¬ 
tration  did  not  reach  C/ C„  of  0.35.  Slow  release  of 
TNT  was  also  shown  by  the  tailing  of  the  desorp¬ 
tion  (right)  side  of  the  breakthrough  results,  Total 
recovery  of  applied  TNT  was  not  obtained .  In  fact, 
on  the  basis  of  the  area  under  the  curve,  we  esti¬ 
mated  that  only  some  40%  of  the  applied  TNT 
pulse  was  recovered  in  the  effluent  solution.  Fig¬ 
ure  9a  illustrates  that,  in  the  presence  of  0,005-M 
Ca(N03)2  as  the  background  solution,  the  mobil- 
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a.  Solid  and  dashed  curves  are  predictions  using  the  transport  equation  with  linear,  Freundlich 
and  Langmuir  retention  models . 


b.  Dashed  and  solid  curves  are  predictions  using  the  transport  equation  with  a  modified 
Langmuir  model  and  a  multireaction  model  (nonlinear  model  with  irreversible  reaction). 


Figure  9.  TNT  breakthrough  results  (closed  circles) for  a  bentonite-sand  column  with  an  input 
TNT  concentration  (C0)  of  10.65  fJg/inL  in  a  0,005-M  Ca(ND^2  background  solution. 
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ity  of  TNT  in  the  bentonite-sand  column  was 
much  more  retarded  than  when  methanol  was 
present  (Fig,  8). 

We  found  that  the  use  of  linear,  Freundlich  or 
Langmuir  retention  with  the  convective-disper¬ 
sion  equation  poorly  described  the  transport  data 
of  Figure  9.  For  all  three  predictions,  the  extent  of 
retardation  or  the  arrival  time  of  TNT  in  the  efflu¬ 
ent  solution  was  overpredicted,  In  addition  to 
delayed  arrival  of  TNT,  the  maximum  predicted 
concentration  was  higher  than  that  experimen¬ 
tally  measured.  The  retention  parameters  associ¬ 
ated  with  the  equilibrium  models  were  based  on 
batch  results  on  the  bentonite  day  (see  Table  3  and 
Fig,  1-3),  It  should  be  emphasized  that  the  three 
models  considered  here  are  equilibrium-type 
models  and  do  not  account  for  irreversible  or 
slowly  reversible  retention  reactions.  As  a  result, 
the  failure  of  these  fully  reversible  equilibrium 
models  to  describe  this  data  set  was  not  surprising, 
In  addition,  the  modified  Langmuir  model  eq  14 
failed  to  describe  the  data  (Fig,  9b), 

Model  predictions  shown  in  Figure  9  illustrate 
the  inadequacy  of  the  equilibrium  models  in  de¬ 
scribing  the  TNT  results  in  a  clay-sand  column  in 
the  presence  of  0,005-M  Ca(N03)2  as  the  back¬ 
ground  solution.  Despite  equilibrium  models  of 
the  Freundlich  and  Langmuir  type  providing  good 
predictions  of  the  batch  data,  one  should  not  auto¬ 
matically  assume  the  validity  of  such  models  to 
transport  experiments.  We  suggest  that  the  batch 
data  of  TNT  on  bentonite  clay  did  not  account  for 
possible  retention  mechanisms  of  the  slow  ad¬ 
sorption-desorption  or  transformation  during  the 
2-day  reaction  time  of  batch  equilibration, 

We  attempted  further  prediction  of  the  TNT 
data  set  of  Figure  9b  using  the  Multi-Reaction  and 
Transport  Model  (MRTM)  described  earlier,  where 
two  mechanisms  were  assumed  to  govern  the 
retention  reactions  of  TNT.  The  first  mechanism 
was  a  fast  reaction  of  the  nonlinear  (Freundlich) 
equilibrium  type  (see  eq  2)  and  the  second  an 
irreversible  mechanism  of  the  first-order  type  (see 
eq  7).  Because  of  lack  of  complete  TNT  recovery  in 
the  effluent  solution  from  the  bentonite-sand  mis¬ 
cible  displacement  column,  incorporation  of  an 
irreversible  reaction  was  necessary  to  account  for 
TNT  phases  that  were  not  In  soil  solution.  Incom¬ 
plete  TNT  recovery  Is  ascribable  to  possible  trans¬ 
formation  reactions  as  well  as  slow  release  of 
sorbed  TNT  phases  on  the  clay  surfaces.  Other 
mechanisms  include  inter-particulate  diffusion  In 
clay  lattices  or  other  inaccessible  sites.  A  diffusion 
mechanism  involving  clays  has  been  postulated 


for  other  ions  by  Sparks  and  co-workers  (see  Sparks 

1989). 

Model  parameters  that  provided  best-fit  pre¬ 
dictions  of  the  TNT  results  of  Figure  9b  are  given 
in  Table  4.  The  MRTM  model  was  capable  of 
predicting  the  peak  concentration  and  the  extent 
of  TNT  retardation  in  the  soil  column.  The  good¬ 
ness  of  fit  of  the  MRTM  model  as  given  by  the 
value  of  r2  was  0.836.  In  addition,  the  standard 
errors  of  model  parameters  were  small,  which 
adds  credence  to  the  applicability  of  the  mecha¬ 
nisms  considered  in  describing  the  Breakthrough 
Curves  (BTCs).  Based  on  model  calculations,  the 
amount  of  TNT  recovery  was  41.9%  following  35 
pore  volumes,  which  compares  well  with  experi¬ 
mental  data.  Estimates  for  Kj  and  b  were  much 
smaller  than  those  based  on  batch  results,  how¬ 
ever. 

Results  for  RDX  concentration  in  the  effluent, 
as  shown  by  the  breakthrough  results  in  Figure  10, 
reveal  that  there  is  a  limited  retention  of  RDX  by 
the  bentonite-sand  column,  This  is  attributable  to 
the  retardation  of  the  BCT,  with  the  arrival  of  RDX 
in  the  effluent  after  two  pore  volumes  and  a  rela¬ 
tive  concentration  (C/Cn)  of  0.5  corresponding  to 
approximately  three  pore  volumes.  In  addition, 
RDX  in  the  effluent  indicates  that  there  is  a  sharp 
rise  of  the  sorption  front  and  on  the  desorption 
side,  which  is  commonly  known  as  self-sharpen¬ 
ing.  This  is  a  strong  indication  that,  for  RDX,  the 
local  equilibrium  is  the  dominant  cause  of  reten¬ 
tion  of  RDX  in  the  bentonite-clay  column, 

The  use  of  an  equilibrium  Freundlich  retention 
with  a  fully  reversible  reaction  in  the  MRTM  to 
describe  RDX  results  is  Indicated  by  the  solid  and 
dashed  curves  shown  in  Figure  10.  Although  the 
predictions  did  not  fully  describe  the  shape  of  the 
curve  for  RDX  in  the  effluent,  they  predicted  the 
arrival  time  of  the  RDX  in  theeffluent .  For  the  solid 
curve,  RDX  predictions  were  obtained  using  non¬ 
linear  least-square  optimization  of  MRTM,  where 
kH  was  set  to  zero.  This  is  because  it  appeared  that 
a  total  recovery  of  the  applied  RDX  pulse  input 
was  obtained  in  the  effluent  solution.  Estimates  for 
the  Freundlich  parameters  Kj  and  b  that  provided 
the  best-fit  of  RDX  results  were  0.8203  cm3/g  and 
0,8048  respectively.  The  errors  of  estimation  for 
both  parameters  were  relatively  small  (see  Table 
4),  It  is  possible  that  the  deviation  from  the  sharp 
rise  of  the  BTC  may  be  caused  by  the  effect  of  a 
large  dispersion  coefficient  or  nonlinearity  of  the 
Freundlich  retention  mechanisms,  It  is  known  that 
a  large  D  value  (or  small  Peclet  number)  or  a  small 
b  value,  or  both,  results  In  a  sharpening  of  the  BTC 
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Table  4.  Best-fit  model  parameters  for  the  TNT  and  RDX  miscible  displacement  experiments. 


Experiment  K,/  k, 


Soil 

number 

Model 

(cm-Vg) 

SE* 

b 

SE 

(hr1) 

SE 

r! 

TNT  pulse 

Bentonite-sand 

102 

LIN 

0,204 

0.009 

— 

— 

— 

— 

0.987 

Bentonlte-sand 

110 

MRTM* 

4,374 

0.073 

0,745 

0.008 

0.072 

0.003 

0.836 

Norwood 

103 

MRTM* 

0.357 

0.019 

— 

0.095 

0,007 

0.978 

Norwood 

105 

MRTM* 

0,400 

0.010 

— 

— 

0.158 

0.003 

0.970 

RDX  pulse 

Bentonite-sand 

112 

FRN 

0,802 

0,036 

0,805 

0,029 

— 

— 

0.919 

Kolin 

112 

MRTM* 

0,725 

0.024 

0,821 

0.014 

0.002 

0,002 

0.971 

+  SE  =  Standard  error. 

*  LIN  ■  linear  model,  MRTM  »  multireaction  and  transport  model,  HRN  =  Preundlleh  model. 


(Selim  et  al.  1990),  We  did  not  attempt  to  adjust  b 
or  D  values  in  this  study,  however,  In  fact,  the  use 
of  a  smaller  b  value  to  produce  a  sharper  front 
requires  further  adjustments  to  the  Kj  value.  Since 
the  D  value  used  is  based  on  that  determined  using 
tritium  breakthrough  results,  we  made  no  attempts 
to  adjust  the  D  value  used  In  the  predictions. 

To  study  the  extent  of  Irreversible  retention  of 
RDX,  we  examined  a  range  of  values  of  the  irre¬ 
versible  rate  coefficient  (kH)  for  transport  predic¬ 


tions  based  on  eq  1 .  Selected  predictions  are  shown 
by  the  dashed  curves  of  Figure  10  for  values  of 
0,01  and  0,02  hr1,  which  resulted  in  10.5  and  21% 
irreversible  retentions  of  RDX.  More  significantly, 
it  is  apparent  that  the  use  of  of  0  produced  the 
best  predictions,  This  is  also  supported  by  a  visual 
estimation  of  the  area  under  the  curve.  Therefore, 
we  conclude  that  RDX  was  highly  mobile  with 
little  retardation  and  that  kinetic  or  slow  retention 
processes  were  absent,  Moreover,  we  did  not  de- 


Pore  Volume  (V/V0) 

Figure  10,  RDX  breakthrough  results  (closed  circles) fora  beutonite-sand  column  until  input 
RDX  concentration  (CJ  of  10,65  Mg/mi  in  a  0,005-M  Cu(N03)2  background  solution,  The 
curves  are  calculations  using  the  transport  equation  with  nonlinear  equilibrium  retention  and 
an  irreversible  reaction  for  ks  o/O, 0,01  and  0.02  hrl 
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Figure  1 1,  RDX  breakthrough  results  (closed  circles)  for  Kolin  uncontamimted  soil  column 
with  input  RDX  concentration  (CJ  of  10,04  pg/tnl  in  0,005-M  Ca(N03)2.  The  solid  curve  is 
calculated  using  the  transport  equation  with  nonlinear  retention  whereas  the  dashed  curve 
includes  irreversible  reaction. 


tect  irreversible  retention  of  RDX  on  bentonite  clay 
surfaces  or  transformation  reactions.  Other  inves¬ 
tigators  (Leggett  1985,  Pennington  and  Patrick 
1990)  have  obtained  low  retention  of  RDX  by  clays 
and  several  soils  during  adsorption-desorption, 
We  investigated  the  mobility  of  RDX  and  TNT 
in  uncontaminated  soil  samples  obtained  from 
near  the  contaminated  site  at  the  Louisiana  Army 
Ammunition  Plant.  Here,  we  introduced  a  p  ulse  of 
a  mixture  of  TNT  and  RDX,  having  concentrations 
of  10,65  and  10.04  ing/L  in  0.005-M  Ca(N03)2  as 
the  background  solution,  to  a  packed  column  of 
the  uncontaminuted  soil  after  a  steady  water  flux 
condition  was  achieved.  We  introduced  approxi¬ 
mately  five  pore  volumes  of  the  TNT  and  RDX 
mixture  to  the  packed  column,  Breakthrough  re¬ 
sults  for  RDX  concentration  in  the  effluent  are 
shown  in  Figure  11  and  indicate  the  arrival  of  RDX 
in  the  effluent  after  about  two  pore  volumes,  where 
peak  (maximum)  RDX  concentration  in  the  efflu¬ 
ent  reached  that  of  the  input  pulse  concentration 
(C/Cu  =  1).  This  was  followed  by  a  gradual  RDX 
concentration  decrease  and  is  indicated  in  the 
figure  by  the  tailing  of  the  desorption  side  of  the 
breakthrough  results,  Comparison  of  data  for  the 
uncontaminated  soil  and  those  for  the  bentonite- 


sand  column  (Fig.  10)  illustrates  the  similarity  of 
RDX  mobility  in  the  two  systems.  Major  differ¬ 
ences  were  found  in  the  overall  shape  of  the  break¬ 
through  results  rather  than  in  the  arrival  time  of 
the  RDX  fronts.  Specifically,  the  RDX  effluent  and 
desorption  sides  for  the  uncontaminated  soil  did 
not  exhibit  the  sharp  rise  (self-sharpening)  seen  in 
the  bentonite-sand  column.  This  may  be  ascrib- 
able  to  a  larger  dispersion  coefficient  (D)  for  the 
uncontaminated  soil  or  the  presence  of  nonlinear 
or  kinetic-type  retention  reactions.  The  extent  of 
RDX  retention  was  similar  for  both  systems,  how¬ 
ever. 

The  use  of  a  nonlinear  retention  of  the  Freundllch 
type  with  no  irreversible  retention  in  the  transport 
model  produced  good  prediction  of  RDX  move¬ 
ment  in  the  uncontaminated  soil,  as  Indicated  by 
the  solid  curve  of  Figure  11.  Parameter  estimates 
for  and  b  that  provided  the  best-fit  of  RDX 
results  were  0,7249 cm3 / g  and  0.8209  respectively. 
The  corresponding  estimation  errors  for  Kj  and  b 
were  small.  When  the  model  was  extended  to 
include  irreversible  retention,  the  estimated  kH  value 
was  0,00212  ±  0.00214  hr1,  which  is  a  poor  esti¬ 
mate.  Moreover,  the  area  under  the  curve  indi¬ 
cated  a  3%  irreversible  retention  only,  Therefore, 


the  RDX  transport  data  shown  for  the  uncontami¬ 
nated  soil  are  consistent  with  those  for  the  bento¬ 
nite-sand  column,  indicating  extremely  high  po¬ 
tential  mobility  of  RDX  in  soils  with  little  or  lim¬ 
ited  retardation. 

Analysis  for  TNT,  using  our  HPLC  procedure, 
did  not  detect  TNT  in  the  effluent  solution  until 
some  six  pore  volumes  eluted  through  the  uncon¬ 
taminated  soil  column  (da  ta  not  shown).  The  TNT 
concentration  did  not  exceed  0,05  mg/L  (or  C/C0 
=  0.005).  After  12  pore  volumes,  the  TNT  de¬ 
creased  and  w  as  no  longer  detectable  by  HPLC,  To 
investigate  the  extent  of  retention  of  TNT  in  this 
soil,  we  continued  application  of  the  background 
solution  (0.005-M  Ca(N03)2)  for  24  pore  volumes. 
Since  no  appreciable  amounts  of  TNT  were  recov¬ 
ered,  we  applied  a  pulse  containing  50%  methanol 
in  the  0.005-M  Ca(N03)2  background  solution  at  a 
constant  flux,  finding  that  TNT  concentrations 
reached  a  maximum  of  0.34  mg/L  at  V/Va  of  2, 
We  observed  no  higher  TNT  concentrations  in  the 
effluent  solution  after  10  pore  volumes  of  50% 
methanol  were  applied;  however,  HPLC  analyses 
indicated  the  presence  of  other  possible  degrada¬ 
tion  products  of  TNT.  Specifically,  a  definable 
peak  on  the  chromatographic  chart  was  consis¬ 
tently  observed  following  the  TNT  peak.  It  Is  pos¬ 
sible  that  TNT  degradation  along  with  strong  re¬ 
tention  reactions  are  responsible  for  the  lack  of 
TNTmobility  in  the  original,  uncontaminated  soil, 

To  further  examine  the  leaching  potential  of 
TNT  in  other  soils,  we  selected  the  Ap  horizon 
from  the  Red  River  valley,  which  is  adjacent  to 
Louisiana  Army  Ammunition  Plant.  The  soil  cho¬ 
sen  was  a  Norwood  soil  (fine  silty,  mixed,  calcar¬ 
eous,  Thermic  Typic  Udlfluvent),  with  79%  sand, 
18%  silt  and  3%  clay.  Breakthrough  results  for 
TNT  (in  0.005-M  Ca(N03)2  background  solution) 
in  two  different  columns  of  Norwood  soil  are 
shown  In  Figure  12.  The  fluxes  of  applied  TNT 
pulses  along  with  soil  moisture  contents  and  bulk 
densities  in  the  soil  columns  were  somewhat  simi¬ 
lar  for  both  breakthrough  results  (see  Table  1 ),  The 
primary  difference  between  the  data  sets  shown  i  n 
the  figure  is  the  concentration  of  TNT  in  the  input 
solution  (C0).  The  shapes  of  the  BTCs  were  gener¬ 
ally  similar.  Specifically,  the  arrivals  of  TNT  In  the 
column  effluent  solutions  were  similar  regardless 
of  C0.  In  addition,  the  shapes  of  the  desorption 
(right)  sides  were  also  similar,  with  a  gradual 
decrease  in  concentration.  The  two  data  sets  also 
indicate  that  TNT  retention  behavior  Is  likely  to  be 
nonlinear  or  kinetic,  or  both,  because  peak  concen¬ 
trations,  as  well  as  the  extent  of  TNT  recovery  in 


the  effluent,  depended  on  Ca, 

The  areas  under  the  BTCs  show  that  less  than 
half  of  the  TNT  applied  into  the  Norwood  soil 
column  was  recovered  when  C0  was  10.64  mg/L 
in  comparison  with  over  60%  recovery  for  C0  of 
100mg/L(seeFig.  12).  Irreversible  retention,  trans¬ 
formation  reactions  or  slow  release  of  TNT,  or  all 
three,  are  possible  mechanisms  that  may  explain 
the  amounts  of  TNT  not  recovered  in  the  effluent. 
Therefore,  we  applied  the  MRTM  to  describe  the 
BTC  data  for  TNT  in  the  presence  of  an  irreversible 
reaction  or  a  sink  term,  along  with  a  fast  (i.e., 
equilibrium)  type  reaction. 

Model  predictions  are  indicated  by  the  solid 
curves  in  Figure  12,  For  both  cases,  we  found  that 
a  linear  approach  (eq  15),  along  with  a  first-order 
irreversible  reaction,  adequately  described  the 
BTCs,  Attempts  to  describe  the  BTCs  using  nonlin¬ 
ear  (Freundlich)  retention  resulted  in  b  values 
close  to  unity  with  relatively  small  standard  er¬ 
rors,  As  a  result,  we  used  the  nonlinear  least- 
squares  optimization  technique  to  fit  the  BTCs, 
where  Xj  and  kH  were  the  only  two  parameters  to 
be  estimated.  The  estimated  Kj  parameter  was 
somewhat  similar  for  both  columns  regardless  of 
C„.  In  contrast,  a  higher  k„  value  was  estimated  for 
the  case  with  the  low  TNT  input  concentration  (C0 
=  10,28  mg/L), 

The  dependency  of  the  irreversible  rate  coeffi¬ 
cient  on  C0  has  been  observed  for  other  solute 
species  by  other  investigators  ( Amacher  et  al ,  1988) , 
Moreover,  for  all  parameters,  the  standard  errors 
of  estimates  were  small.  Based  on  model  calcula¬ 
tions,  the  amounts  of  TNT  in  the  slowly  reversible 
or  irreversible  phase  were  estimated  as  56.5  and 
33.8%  of  that  applied  for  C0  of  10.28  and  100 
mg/L  respectively. 

The  model  results  in  Figure  12  indicate  that  it 
was  incapable  of  describing  the  slow  release  of 
TNT  during  leaching,  regardless  of  the  TNT 
concentration  in  the  pulse  input  C0.  Kinetic  bchav- 
Iorof  a  chemical,  biological  or  physical  nature  may 
be  responsible  for  the  slow  release  of  TNT  in  the 
effluent,  Moreover,  the  reason  for  such  high  reten¬ 
tion  capacity  of  TNT  in  Norwood  soil  is  not  com¬ 
pletely  understood.  Norwood  soil  has  a  high  pH 
with  low  cation  exchange  capacity,  a  low  clay 
content  and  an  extremely  low  organic  matter  con¬ 
tent,  We  thought  that  TNT  would  behave  as  a 
nonreactive  soiute,  with  little  if  any  TNT  retained 
irreversibly  within  the  soil  column.  Therefore, 
additional  investigations  are  needed  to  determine 
the  behavior  of  TNT  in  soils  such  as  the  Norwood 
soil  under  flow  conditions. 
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Relative  Concentration  (CJCJ  Relative  Concentration  (C/C^ 


Pore  Volume  (V/V0) 

a .  Input  concentration  (CJ  oflOOpg/mL  in  0,005-M  Ca(NO 3)2. 


b,  input  concentration  (C0)  of  10,28  pg/mL  in  0.005-M  Ca(NO^)2. 

Figure  1 2.  TNT  breakthrough  results  fora  Norwood  sail  column.  The  solid  curve  is  the  prediction 
using  the  transport  equation  with  a  linear  retention  and  irreversible  reaction. 


19 


Table  5,  Concentrations  (jag/g)  of  explosive  residues  released  from 
contaminated  (Kolin)  soil  samples  from  the  Louisiana  Army  Am¬ 
munition  Plant. 


AmDNT 

HMX  RDX  TNII  TNT  sample 


Sample  A  39.8  175  18,4  164  30,6 

Sample  B  36,3  127  32,4  76,1  35,2 


Contaminated  soil 

We  studied  the  release  of  TNT  and  RDX  from 
the  contaminated  (Kolin)  soil  samples  (A  and  B) 
obtained  from  the  Louisiana  Army  Ammunition 
Plant.  The  uncontaminated  (original)  samples, 
along  with  duplicate  uncontaminated  samples, 
were  analyzed  by  CRREL  scientists  for  the  total 
contaminants, 

The  results  for  HMX,  RDX,  TNB,  TNT  and 
AmDNT  are  shown  in  Table  5,  Contaminated  soil 
A  contained  considerably  higher  amounts  of  RDX 
and  TNT  than  sample  B,  The  reasons  for  such 
differences  are  not  known  and  may  be  a  reflection 
of  the  spatial  heterogeneity  of  contaminant  con¬ 
centration  in  the  site.  Uncontaminated  samples 
did  not  show  detectable  amounts  of  TNT  or  RDX. 

Results  of  the  re’ease  study  from  the  contami¬ 
nated  samples  are  3hown  in  Figure  13.  The  RDX 
and  TNT  release  data  represent  averages  of  three 
replications.  For  each  extraction  step,  i.e.,  follow¬ 
ingdecanting  and  additions  of  CafNC^fo  solution, 
the  time  of  release  reaction  was  1  day.  For  sample 
A,  the  release  of  RDX  and  TNT  continued  for  8 
days  until  concentrations  below  the  detection  lim¬ 
its  of  HPLC  measurements  were  reached.  Integra¬ 
tion  of  the  total  amount  released  Indicated  that 
more  than  90%  of  the  TNT  and  RDX  release  took 
place  within  the  first  3  days  of  reaction.  This  obser¬ 
vation  was  consistent  for  both  samples  A  and  B. 
Differences  in  the  rates  of  release  are  perhaps 
caused  by  the  total  initial  TNT  and  RDX  being 
considerably  higher  in  sample  A.  The  concentra¬ 
tion  levels  were  lowered  by  two  orders  of  magni¬ 
tudes  after  4  to  5  days  of  continued  release-extrac¬ 
tion. 

An  important  aspect  of  release  studies  Is  that 
one  can  estimate  a  mass  balance  of  the  contamina¬ 
tion  at  different  release  times.  Based  on  our  release 
results,  the  total  TNT  amounts  released  were  103 
and  42  pg/g  soil  for  samples  A  and  B  respectively 
(on  the  basis  of  an  average  of  the  three  replications 
of  concentration  at  each  release  step),  These 
amounts  represent  60%  of  the  total  initial  amounts 


present  in  the  contaminated  samples.  On  the  other 
hand,  for  RDX,  the  amounts  released  were  215  and 
155  pg/g  soil  for  samples  A  and  B  respectively — 
amounts  that  were  somewhat  higher  than  the 
Initial  values. 

Transport  of  TNT  and  RDX  in  contaminated 
soil  samples  A  and  B  is  shown  in  Figure  14.  To 
closely  represent  the  ionic  strength  of  the  soil 
solution  in  the  top  soil  layers,  the  background 
solution  applied  was  0.005-M  Ca(N03)2,  The  re¬ 
sults  clearly  show  an  initial  rise  in  concentrations 
of  both  TNT  and  RDX  in  the  effluent  solution 
because  the  contaminated  soils  were  initially  air- 
dried  ,  As  the  background  solution  was  introduced 
to  the  dry  soil  and  water  movement  in  the  unsat¬ 
urated  state  commenced,  the  solutes  present  (TNT 
and  RDX)  dissolved  in  the  soil  solution  as  the 
wetting  front  advanced  in  the  soil  column.  This 
process  continued  until  the  wetting  front  reached 
the  end  of  the  column  and  effluent  was  collected. 

The  results  of  Figure  14  show  higher  concentra¬ 
tions  of  TNT  and  to  a  lesser  extent  RDX  from 
contaminated  soil  sample  A  than  from  sample  B. 
These  results  are  consistent  with  the  release  data 
from  the  batch  studies  and  reflect  total  higher 
initial  amounts  of  TNT  in  contaminated  soil  A. 
From  the  area  under  each  BTC,  we  also  found  that 
100%  of  the  RDX  initially  present  In  each  soil 
sample  was  released.  In  fact,  only  10  pore  volumes 
of  applied  background  solution  were  required  for 
approximately  90%  of  the  RDX  contamination  to 
leach  out  from  the  contaminated  soil  column,  These 
results  are  also  consistent  with  earlier  findings  for 
the  bentonite-sand  and  Norwood  soil  columns, 
indicating  that  RDX  Is  not  strongly  retained,  As  a 
result,  we  conclude  that  RDX  Is  strongly  suscep¬ 
tible  to  leaching  losses  to  lower  soil  depths. 

Release  and  movement  of  TNT  in  the  contami¬ 
nated  soil  columns  Indicate  a  slow  release  process 
and  a  gradual  decrease  in  effluent  concentration 
(Fig.  13).  Such  a  slow  release  behavior  is  often 
associated  with  kinetic  retention  of  the  linear  and 
nonlinear  mechanisms,  In  contrast  to  RDX,  the 
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b.  Sample  B, 

Figure  13.  Desorption  or  release  of  RDX  and  TNT  from  the  contaminated  soil  samples 
obtained  from  the  Louisiana  Army  Ammunition  Plant.  Each  extraction  step  represents  1  day 
of  desorption-release  time. 
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b.  Sample  B. 


Figure  14,  TNT  and  RDX  breakthrough  results  (pore  volume,  V/V0,  vs.  relative  concentration, 
C/C0)  for  the  contaminated  soil  samples  obtained  from  the  Louisiana  Amy  Ammunition  Plant. 
The  input  solution  was  0.005-M  Ca(N03)2. 
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extent  of  TNT  recovery  was  considerably  less  than 
100%,  In  fact,  after  some  40  pore  volumes,  only 
50%  of  TNT  in  contaminated  soil  A  and  65%  of  that 
in  soil  B  were  leached  out. 


SUMMARY  AND  CONCLUSIONS 

The  objective  of  this  investigation  was  to  deter¬ 
mine  the  potential  mobility  of  explosive  contami¬ 
nants  in  the  soil  to  the  groundwater.  To  achieve 
this,  TNT  and  RDX  sorption-desorption  and  trans¬ 
port  experiments  were  conducted  on  a  reference 
clay,  contaminated  soil  from  the  Louisiana  Army 
Ammunition  Plant  and  two  uncontaminated  soils, 
Standard  batch  experiments  were  used  to  deter¬ 
mine  the  extent  of  adsorption  and  release  of  TNT 
for  a  wide  range  of  concentrations.  In  addition, 
miscible-displacement  techniques  were  used  to 
study  the  transport  of  TNT  and  RDX  in  soil  col¬ 
umns  under  steady  water  flow.  We  found  that 
TNT  isotherms  for  bentonite  clay  after  1  day  of 
sorption  can  be  described  equally  well  with  linear, 
Freundlich,  Langmuir  and  modified  Langmuir 
models.  Results  of  several  desorption  steps  from 
successive  dilutions  Indicated  that  both  adsorp¬ 
tion  and  desorption  isotherms  coincided  closely. 
This  indicated  a  lack  of  hysteretic  behavior  of  TNT 
in  bentonite  clay.  Mass  balance  calculations  indi¬ 
cated  that  TNT  retention  is  also  fully  reversible.  As 
a  result,  for  the  bentonite  clay  used  in  this  study, 
there  were  no  irreversible  interactions  of  TNT, 
such  as  degradations,  within  the  time  of  the  batch 
experiments. 

Transport  results  from  bentonite-sand  miscible- 
displacement  columns  indicated  that  TNT  was  a 
highly  mobile  contaminant  and  was  considered 
fully  conservative  in  the  presence  of  methanol  as 
the  background  solution.  In  fact,  TNT  transport 
closely  resembled  that  of  a  nonreactive  spei  ies  in 
porous  systems.  In  contrast,  TNT  was  strongly 
retained  (more  than  50%  of  the  amount  applied) 
within  the  bentonite-  sand  column  for  a  prolonged 
period  when  the  applied  TNT  pulse  was  in  0.005- 
M  Ca(N03)2  as  the  background  solution.  The  use 
of  a  cc  tvection-dispersion  equation  with 
Freundlich  (equilibrium)  retention  and  a  first-or¬ 
der  irreversible  reaction  adequately  predicted  the 
TNT  transport  results,  Poor  model  predictions 
were  obtained  when  best-fit  parameter  values, 
based  on  batch  retention  results,  were  used  to 
describe  the  behavior  of  TNT  during  transport  in 
the  bentonite-sand  column.  Reasons  for  this  fail¬ 
ure  were  not  completely  understood  and  may  be 


attributable  to  kinetic  or  time-dependent  TNT  re¬ 
actions  as  well  as  physical  retention  by  diffusion 
into  clay  lattices. 

Transport  results  from  Norwood  (fine  silty) 
soil  columns  revealed  that  TNT  was  strongly  re¬ 
tained  in  this  low-organic-matter  (0.32%)  and  low- 
clay-content  (3%)  soil  in  a  manner  similar  to  tha  t  in 
the  bentonite-sand  columns.  We  found  that  the 
use  of  a  linear  reversible  retention  along  with  a 
first-order  irreversible  (sink)  provided  a  good  de¬ 
scription  of  breakthrough  results  regardless  of 
TNT  input  pulse  concentrations  (C0's).  However, 
the  observed  slow  release  of  TNT  from  the  mis¬ 
cible  displacement  columns  was  not  well  pre¬ 
dicted  by  the  model,  nor  was  it  observed  for  the 
bentonite-sand  column,  Slow  release  is  indicative 
perhaps  of  kinetic  retention  of  TNT  in  this  soil. 

RDX  and  TNT  were  released  from  the  contami¬ 
nated  Kolin  soil  for  8  days  before  reaching  concen¬ 
trations  below  the  detection  limits  of  experimental 
(HPLC)  measurements,  In  excess  of  90%  of  the 
TNT  and  RDX  from  the  contaminated  soil  was 
released  within  the  first  3  days  of  reaction.  Based 
on  our  results  from  two  subsamples  and  three 
replications,  the  amount  of  TNT  released  from  the 
contaminated  soil  represented  60%  of  the  amount 
originally  present.  In  contrast,  RDX  was  com¬ 
pletely  released.  This  finding  supports  transport 
data  from  the  miscible-displacementexperiments, 
which  indicated  that  100%  of  RDX  in  the  contami  ¬ 
nated  sample  was  released  when  a  solution  of 
0,005-M  Ca(N03)2  was  continuously  applied  for 
20  pore  volumes.  The  extent  of  TNT  recovery  was 
considerably  less  than  that  for  RDX.  In  fact,  after 
some  40  pore  volumes  of  continuous  leaching, 
about  50%  of  TNT  remained  in  the  contaminated 
soil  samples. 
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